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Abstract 

Kinetically constrained models (KCMs) are models of glass formers based on the concept of dynamic 
facilitation. This concept accounts for many of the characteristics of the glass transition. KCMs usu- 
ally display a combination of simple thermodynamics and complex glassy dynamics, the latter being a 
consequence of kinetic constraints. Here we show that KCMs can be regarded as systems whose config- 
uration space is endowed with a simple energy surface but a complicated geometry. This geometry is 
determined solely by the kinetic constraints, and determines the dynamics of the system. It does not 
affect the overall distribution of states. Low temperature dynamical slow-down is then a consequence of 
the competition between the paths allowed by the geometry of configuration space, and those leading to 
energy relaxation. This competition gives rise to dynamical bottlenecks unrelated to static properties. 
This view of the dynamics is distinct from that based on an underlying static rugged energy landscape. 
We illustrate our ideas with simple examples. 

1 Introduction 

The view of the glass transition (see |2J |3j for general reviews) as a dynamical phenomenon emerges 
naturally from the real-space picture of supercooled liquids developed in 01|S1|SJ[7||S]. This work demonstrates 
that many of the properties of supercooled liquids can be simply explained by the existence and statistical 
properties of spatially correlated dynamics, or dynamical heterogeneity pi llUlITT] . Dynamical heterogeneity 
has been observed both in deeply supercooled liquids 12 , in the mildly supercooled regime in simulations 
[T5| . and in experiments on colloidal suspensions [T%| . 

Heterogeneous dynamics, in combination with structural homogeneity, as observed in glass formers, 
appears in systems possessing both localized excitations and facilitated dynamics llfil H7j . The concept 
of dynamic facilitation stipulates that the dynamics of a local region is suppressed if is bordered by immobile 
regions, and thus reflects the idea that particles in a jammed liquid cannot rearrange unless adjacent to 
mobile particles. Dynamic facilitation is realized in KCMs, which are simple models of non interacting 
particles or excitations in a lattice whose dynamics is subject to kinetic constraints (see |18j for a review). 
Examples of KCMs are spin facilitated models, like the Fredrickson-Andersen (FA) model ^7] or the East 
model ^!5]> an( i constrained lattice gases, such as the Kob- Andersen model [2TJ| . 

In this paper we study the effect that dynamic facilitation has on the geometry of configuration space, 
and how in turn this affects the dynamics |23- We show, by generalizing KCMs to continuous degrees of 
freedom, that kinetic constraints impose a metric structure on configuration space which does not affect the 
static distribution of states, and therefore the thermodynamics, but which has a dramatic effect on the paths 
between states at low energies. We show that dynamical slow-down at low temperatures occurs because the 
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system is forced to evolve along configuration space paths of least metric distance, or geodesies, which are 
distinct from the paths that relax the energy. This competition between available trajectories and energy 
relaxation gives rise to dynamical barriers, or bottlenecks, which are unrelated to the statics of the system. 
This view of the dynamics is distinct from that based on the idea that a static rugged energy landscape is 
responsible for the properties of glass formers (see e.g. [3]). 

This paper is organised as follows. In section |3 we discuss the geometric interpretation of dynamic 
facilitation. In sections and 0] we illustrate our ideas with simple examples. In section we state our 
conclusions. 



2 Geometrical interpretation of constrained dynamics 

The real-space picture of supercooled liquids of Refs. El is based on two observations. First, at 
low temperature very few particles in the liquid are mobile, and these mobility excitations are localized in 
space. Second, regions of the liquid cannot become mobile unless their neighbouring regions are mobile. 
These observations are implemented as follows. A supercooled fluid in d spatial dimensions is coarse-grained 
into cells of linear size of the order of the static correlation length as given by the pair correlation function. 
Cells are classified by a scalar mobility field, x, identified by coarse-graining the system on a microscopic 
time scale. Mobile regions carry a free energy cost, and when mobility is low interactions between cells 
are not important. Adopting a thermal language, we expect static equilibrium to be determined by the 
non-interacting Hamiltonian, 

N 

H[x] = Y l M^), (1) 

where /i — 1, . . . , N labels the lattice site, is the mobility field at site /i, and / M is an analytic function of 
x^ chosen so that the ground state of the model is such that | is small. At low mobility, the distinction 
between single and multiple occupancy is irrelevant, as is the distinction between discrete or continuous 
degrees of freedom. In what follows we assume x^ to be real. The dynamics of the mobility field is given by 
a master equation, 

d t P(x,t) = J2 C »( x ) £»P(x,t), (2) 

where P (x, t) is the probability that the system has configuration x = {x^ \ /j, = 1, 2, N} at time t. The 
local operators £ M are the same as those for unconstrained local dynamics. Their action on P can be written 

4, p(x) = j2 w ( x '^ -» p o=' M > {x^}) - J2 w ^ x>1 x '^ p ( x>1 > K 5 ^}) > ( 3 ) 

where w{x^ — > x'^ 1 ) is the probability of going from configuration (x^, {x"^}) to configuration {x"^ 11 }) 
in unit time. We ensure a unique equilibrium configuration exists by requiring (j2J) to obey detailed balance 
with respect to QJ, i.e. 

W ( X ^ X "> = e -p(H[x']-H[x])_ / 4 s 

w(x' — > x) 

The kinetic constraint, C M (x), is designed to suppress the dynamics of cell /j, when surrounded by immobile 
regions. It must also be such that Eq. J2J) satisfies detailed balance, which is achieved if C M (x) does not 
depend on x M itself. To reflect the local nature of dynamic facilitation we allow C M to depend only on the 
{x v } of nearest neighbours of /i, and require that is small when local mobility is scarce. Equations 0-0 
define a generic KCM. 
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We now show that Eq. has a simple geometrical interpretation. Using standard methods [221 we P ass 
from the master equation J5J to the Fokker-Planck equation 

3 t P(M)=£4 FP >PM), (5) 

where 

The prime on denotes differentiation with respect to x' 1 . Equation © has a simple physical interpretation: 
it describes single-particle driven diffusion on an TV-dimcnsional curved space {x^ | /i = 1,2, ...,N} . The 
form of the driving term in (JSJ comes from the Hamiltonian and the kinetic constraint. But the curvature 
of this space is due solely to the kinetic constraint: from the form of a Fokker-Planck equation on a curved 
manifold 22; we identify the inverse metric tensor of this space as g^"(x) = 5^C^ (x). Thus the metric reads 

9nu(x) = S^C- 1 0) . (7) 

No summation is implied. The effect of dynamic facilitation is to endow the configuration space x with a 
nontrivial metric. Since dynamic facilitation is responsible for the interesting dynamical properties of these 
models, we expect that these dynamical properties can be extracted from the metric tensor. 

The effect of the metric on configuration space trajectories is made clearer by writing a formal solution 
to equation © 23j: 

r b 

P{x b ,t b \x a ,t a ) = y/\ji\Vxe- s W, (8) 

J a 

where P(x b , tb\x a , t a ) is the conditional probability that a system starting at point x a in configuration 
space at time t a is found at point x b at time t b . We have introduced the determinant of the metric, 
\g\ = det [gfiu) = Il^/T 1 - The path integral © is taken over all configurations of the system, weighted by 
the dynamic action S[x\. Up to constant terms this action reads 

S[x] = i (H[x b ] - H[x a \) + I dt (x» 9llv x» + ||<r H) , (9) 

where dots denote differentiation with respect to time, and the summation convention for once-repeated 
upper and lower indices is implied. 

The action JSJ) weights the configuration space trajectories of our system. The system's classical tra- 
jectory, i.e. the path taken by the dynamics when T — > 0, is one of extremum S. At finite T, stochastic 
fluctuations cause perturbations about this path. Equation (|5J shows clearly the way the various ingredi- 
ents of a dynamically facilitated model influence its trajectory through configuration space. The first term 
depends only on the Hamiltonian evaluated at the endpoints. It gives rise to Boltzmann weights in 
the path integral (JHJ) and plays no role in choosing the trajectory of the system. The path is chosen by a 
competition between the two terms in the integral of Eq. ©. The first term, S± = T^ 1 J dtx^g^x" , is 
purely geometrical. It is the action for free motion in a curved background. It depends on the metric g^, 
and therefore on C _1 . The second term, 52 = T _1 J dt d^H g^ u d u H , depends on the Hamiltonian (JTJ, but is 
proportional to the inverse metric g^ v , and so to C. At low temperatures, we expect 52 to be important only 
far from equilibrium, when mobility is plentiful — or equivalently, when the magnitudes of the coordinates 
\x^\ are large and C is large — and the system's relaxation will take place by way of gradient descent on the 
energy surface. When mobility has diminished significantly the constraints C become small. We then expect 
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iSi to dominate. Hence the classical trajectory of the system in the region of weak dynamic facilitation will 
be governed by extremum Si . This corresponds to the equation of a geodesic 24 : 

x» + T^x a xP = 0. (10) 

The Christoffel symbols T^a are defined as = \g^ v {g V a,p + 9pv,a — 9af3,v), and the metric tensor raises 
and lowers indices. Equation IjlQI) tells us how the kinetic constraint controls the system's trajectory in 
configuration space in the dynamically constrained regime. When mobility is scarce, the constraint C is 
small, and the metric distance between points in configuration space, s a b = J \J g ^dx^ dx v , is large, by 
virtue of Q. By minimising this distance, and so following Eq. (|1(J|) . we expect the system to follow a 
path that does not lead directly to a configuration of lowest energy. This is a geometrical interpretation of 
dynamical arrest. 



3 A simple example 

We now illustrate these ideas with a simple example. We define a model of two variables, {x^} — {x, y}, with 
each variable constraining the other. Although such a model cannot exhibit true glassiness it is useful because 
we can visualise directly the two-dimensional configuration space it inhabits. We choose the Hamiltonian 
Q to be quadratic, 

H = ±(x 2 +y 2 ), (11) 

where J sets the energy scale. We choose the constraint functions to be quadratic also, with each coordinate 
constraining the other: 

r — n, 2 I 

Thus, the Fokker-Planck equation (JSJ is determined. Equation Q gives the metric of this system as 



C x =y 2 , C y = x 2 . (12) 



^ = (Yi/°0- (13) 

The metric is singular, and so distances in this configuration space become very large, when x or y vanish, re- 
flecting the tendency of the kinetic constraint to suppress the dynamics of the system in the regime of low mo- 
bility. The space described by this metric is negatively curved with curvature scalar R = —4 {x 2 /y 2 + y 2 /x 2 ) 
|24j . The geodesies follow from (|10(1 which in this case read 

y'x + (y/x) 3 y 2 - 2xy = 0, xy - (x/y) 3 x 2 - 2xy = 0. (14) 

To illustrate the effect of the kinetic constraint on the trajectories chosen by the system we first cast Equation 
JSJ) into a form more convenient for numerical simulation. In the usual way |25| we recognize that © is 
equivalent to the coupled Langevin equations 



x{t) \_ f Jy 2 (t)x(t) \ , / \y(t)\ m (t) 



y{t) ) 1 J * ■(*)!/(*) J V N(*)M*) 



(15) 



where the Gaussian white noise variables r)i{t) have zero mean and variance {fii{t)r}j (t 1 )) — 2T6(t — t')5ij. 
In this representation the kinetic constraint manifests itself as a state dependent rate: the driving term is 
linear in C while the noise term goes as y/C. The unconstrained version of the same model, i.e. one for which 
(C x ,C y ) = (1,1), would have Langevin equation x M (t) = ~Jx^ + rj^it). 
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We now show explicitly that low temperature trajectories with fixed initial and final configurations are 
approximated by the geodesic equation 1)14(1 ■ We have used Transition Path Sampling (TPS) J2U to ensure 
that we obtained a set of trajectories which both obeyed the endpoint conditions are were well-sampled from 
the dynamical action 10 ■ In order to generate out-of-equilibrium trajectories for the relaxation from high 
energy configurations we used the Crooks-Chandler |27] algorithm; alternatives such as the local algorithm 
of Ref. |2EJ are also suitable. 

Figure ^ presents the comparison of the dynamical evolution of our toy KCM with the geometry of its 
configuration space. TPS was used on the system defined by Eq. (|15f) to generate trajectories between the 
configurations indicated by a star. In Fig. la we show both a characteristic trajectory (thin red path) and 
the average over the set of all trajectories obtained via TPS (thick blue path). We also show geodesies 
(black curves) joining the final configuration to the initial configuration (indicated by a star), and to other 
configurations along the gradient descent path (indicated by black squares). The solutions to lfH|) are 
degenerate: there is one geodesic going left from the initial point and another going right. The distinction 
depends on the initial conditions on the path. In the TPS simulations we force the initial displacement to 
be to the left to obtain a set of trajectories which follows the left hand geodesies. These are the paths shown 
in Fig. la and lb. In Fig. lb we also show the average over a small set of right hand TPS trajectories and 
the corresponding pairs of geodesies. 

In Fig. ^ we see the features described in the previous section. Initially the system is in the high energy 
region and it relaxes approximately by gradient descent. As it reaches the region of weak dynamic facilitation 
at low energies, the system relaxes by following approximately the geodesies of configuration space. Shown 
for comparison in Fig. la are the trajectories followed by an unconstrained version of the same model (dashed 
green path). In this case relaxation is by gradient descent. Given that the rate for motion in the x (resp. 
y) direction vanishes when y (resp. x) vanishes, trajectories have to cross the coordinate axis orthogonally, 
which gives trajectories (and geodesies) a "rectangular" shape. Equipotential lines are concentric circles, 
and the conflict between allowed dynamical paths and energy relaxation is evident. 

4 An even simpler example 

We can gain more insight into the nature of dynamic facilitation by considering an even simpler, one- variable 
model. This can be thought of as modeling facilitated dynamics for an order parameter, rather than for 
microscopic variables. Equations JHJ and tell us that an TV-variable, dynamically constrained system looks 
like a single particle diffusing in N space dimensions. Thus a one-variable dynamically constrained model 
looks like a single particle diffusing in one space dimension. We can exploit the fact that one dimensional 
spaces are always flat [23] to 'remove' the kinetic constraint. We will see below that this can be done at the 
expense of introducing an effective free energy barrier. 

Consider a model with one degree of freedom, x, with Hamiltonian H = Jx 2 /2. We choose the 
dynamical constraint to be C(x). This differs from the higher-dimensional cases in that the variable constrains 
itself. Consequently, the transformation from Fokker-Planck to Langevin form is ambiguous |25| . and depends 
on whether one choses pre-point or mid-point discretisation of time-dependent quantities (Ito or Stratonovich 
calculus). We choose the former, so that our model is causal. The Langevin description of our one- variable 

model is then 

x(t) = -JxC(x)+TC'(x) + y^C(xjri(t), (16) 

where the prime denotes differentiation with respect to x, and the noise term r](t) has zero mean and 
variance (rj(t)rj(t')) = 2T5(t — t'). Since in this case the constraint depends on the variable it constrains, 
the second term on the right hand side of (|16J) is needed to ensure an equilibrium density proportional to 
exp (— Jx 2 /2T) [2H]. In the higher-dimensional cases we considered above there was no self-facilitation and 
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Figure 1: Configuration space trajectories for the toy model described in the text. Panel (a): T = 2.0, J = 
1, t b -t a = 129, {x a , Va) = (3.6, 4.0), (x b , y b ) = (-0.6 ± 0.01, -0.5 ± 0.01). We used TPS with the Crooks- 
Chandler 27[ algorithm to generate 2.78 x 10 8 trajectories, starting from point a (indicated by a star at the 
top right), of which 11403 landed in the end-point b (indicated by a star at the bottom left). The thin red 
path is a representative trajectory. The thick blue path is the mean of those trajectories. The dashed green 
path is the mean of 273 trajectories reaching b from a for an unconstrained version of an otherwise identical 
model. Numerical solutions to the geodesic equation l|14|) joining b to three other points (shown as black 
squares) are overlaid. The gray circle encloses the equilibrium region (x 2 + y 2 ) = T/J. Panel (b): a similar 
plot starting from (x a ,y a ) — (1.6,2.0). Again, 2.78 x 10 s trajectories were generated. The thick blue path 
shows the mean of 2199 trajectories landing in region b which follow the left hand side route. The average 
path over 6 right hand side trajectories is also shown in green. 
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this term was absent. The corresponding unconstrained model has C{x) = 1, and its Langevin equation is 
x{t) = -d x {Jx 2 /2)+r,{t). 

The appearance of the multiplicative noise term yJC{x)rj{t) in Eq. Ijl6(l is due to the metric-augmented 
diffusion term in JJ]J, and is therefore a consequence of the kinetic constraint. The statement that a one- 
dimensional space is flat is equivalent to the statement that the noise may be made additive by a local change 
of variables. If we make the change of variables 

^/"JIt <17) 

we can use Ito's formula [2EJ to rewrite (|16[l as 

y{t) = -Jx^/C{x) + ^C{x) d x \nC(x) +fj(t), (18) 

where x should be eliminated in favour of y. The noise fj(t) has the same mean and variance as rj(t). Provided 
we can write the deterministic terms in (118(1 as the gradient of an effective potential we can regard Equation 
(|18Jl as unconstrained, in the sense that the prefactor of the noise does not contain the dependent variable. 
Consider the specific case of a quadratic constraint, C(x) — x 2 . Equation (|To)) becomes 

x = -Jx 3 + 2Tx + xrj, (19) 

where we have assumed x > 0. The change of variables (|17(l now amounts to a Cole-Hopf transformation, 
y = Itlx. Hence the regime where dynamic facilitation is relevant corresponds to y <C 0. Equation (|18[) now 
reads 

y = -d y (±e 2 y-Ty)+fj. (20) 

Equation (|2U[) looks like an unconstrained model with an effective entropic barrier AF = —Ty = —T Ins, 
which prevents the particle from accessing the dynamically constrained regime ofy<CO=>a;~0 |3(J| . Figure 
121 illustrates this. 



5 Conclusions 

In this paper we have presented a geometric interpretation of dynamic facilitation. By generalizing kineti- 
cally constrained models to continuous degrees of freedom we have shown that in these systems the kinetic 
constraints can be seen as endowing configuration space with a non-trivial metric structure. This metric 
structure in turn determines the dynamical trajectories. In the regime where the kinetic constraint is domi- 
nant, i.e. at low temperatures or energies, the geometry forces the dynamics to proceed through paths close 
to geodesies, rather than close to paths which relax the energy. This conflict between available paths and 
energy relaxation gives rise to dynamical bottlenecks or barriers, despite the fact that there are no static 
barriers in these systems. 

Our view of dynamical arrest in terms of the geometry of configuration space is distinct from that based 
on the idea of motion in an effective rugged energy (or free energy) landscape (see 3 and references therein). 
In the landscape perspective it is often assumed that relevant features of the energy surface, like local minima 
and saddle points, determine the behaviour of the system, and that knowledge of the statistical properties 
of these features is sufficient to explain both dynamic and thermodynamic properties. In our approach the 
effective energy surface displays no particular features, the thermodynamics is uninteresting, and there are 
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Figure 2: An illustration of how a constrained problem may have an alternative representation as an 
unconstrained one with an extra free energy barrier (from the example of section QJ. The constraint C(x) 
fixes the metric g(x) = C^ 1 (x) of space x. This implies that metric distances become large in the region of 
small x. Under the transformation Eq. I|17[) the problem changes into a non constrained one, g{y) — 1, but 
with an extra free energy barrier AF(y) which excludes the system from the low energy region of large \y\. 

no special configurations, such as local minima or saddles, which play any major roles either in the statics 
or in the dynamics. All the interesting structure is in the paths between configurations, and therefore in 
the dynamics. In this picture it is the effective metric structure of configuration space that distinguishes a 
glassy system from a normal one, and not the ruggedness of its energy surface. 

We are grateful to D. Chandler for discussions. We acknowledge financial support from EPSRC Grants 
No. GR/R83712/01 and GR/S54074/01, the Glasstone Fund, and Linacre College Oxford. 
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